Astronomy & Astrophysics manuscript no. 7320 


© ESO 2008 


February 1, 2008 





Letter to the Editor 



The IC1396N proto-cluster at a scale of -250 AU * 

R. Neri 1 , A. Fuente 2 , C. Ceccarelli 3 , P. Caselli 4 - 5 , D. Johnstone 6 ' 7 , E.F. van Dishoeck 8 , F. Wyrowski 9 , M. Tafalla 2 , B. 

Lefloch 3 , and R. Plume 9 



o 
o 

(N 



oo 



Oh- 

6 



> 

O 

o 



Institut de Radioastronomie Millimetrique, 300 rue de la Piscine, 38406 St Martin d'Heres Cedex, France 

Observatorio Astronomico Nacional (OAN), Apdo. 112, E-28803 Alcala de Henares (Madrid), Spain 

Laboratoire d'Astrophysique de l'Observatoire de Grenoble, BP 53, 38041 Grenoble Cedex 9, France 

Osservatorio Astrofisico di Arcetri (INAF), Largo E. Fermi 5, 50125 Firenze, Italy 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 0213 

Department of Physics and Astronomy, University of Victoria, Victoria, BC V8P 1A1, Canada 

National Research Council of Canada, Herzberg Institute, 5071 West Saanich Road, Victoria, BC V9E 2E7, Canada 

Leiden Observatory, PO Box 9513, 2300 RA Leiden, Netherlands 

Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

University of Calgary, 2500 University Drive NW, Alberta T2N 1N4, Canada 



Received February 17, 2007; accepted April 12, 2007 



ABSTRACT 



Aims. We investigate the mm-morphology of IC 1396 N with unprecedented spatial resolution to analyze its dust and molecular gas 
properties, and draw comparisons with objects of similar mass. 

Methods. We have carried out sensitive observations in the most extended configurations of the IRAM Plateau de Bure interferometer, 
to map the thermal dust emission at 3.3 and 1.3mm, and the emission from the 7=13^ — » 1 2^ hyperfine transitions of methyl cyanide 
(CH 3 CN). 

Results. We unveil the existence of a sub-cluster of hot cores in IC 1396N, distributed in a direction perpendicular to the emanating 
outflow. The cores are embedded in a common envelope of extended and diffuse dust emission. We find striking differences in the 
dust properties of the cores (J3 - 0) and the surrounding envelope (j3 - 1), very likely testifying to differences in the formation and 
processing of dust material. The CH 3 CN emission peaks towards the most massive hot core and is marginally extended in the outflow 
direction. 
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1. Introduction 

Intermediate-mass young stellar objects (IMs) (protostars and 
Herbig Ae/Be stars with M+ ~2-10M o ) are crucial to star for- 
k> \ mation studies because they provide the link between evolu- 
5_j ■ tionary scenarios of low- and high-mass stars. These objects 
& ' share many similarities with high-mass stars (clustering, PDRs). 
However, to study them presents decided advantages, compared 
to massive star forming regions, as many of them are located 
close to the Sun (<1 kpc) and in regions of reduced complexity. 

The lack of spatial resolution is always a potential source of 
confusion in the identification of isolated stellar structures. This 
issue becomes increasingly important in the case of low-mass 
and intermediate-mass stars, which are predisposed to form in 
clusters. While much is now known about isolated star forma- 
tion, the onset of cluster formation remains less clear and sys- 
tematic surveys are needed to gain a deeper understanding (Testi 
et al. 2000). Aperture synthesis in the millimeter range provides 
excellent means to investigate and characterize the earliest stel- 
lar clustering. 
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In this Letter, we present interferometric continuum and 
molecular line observations of the IM protostar IC 1396 N at the 
highest spatial resolution of the PdBI (Karastergiou et al. 2006). 
IC 1396N is a -440 L Q (Sugitani et al. 2000) protostar located 
at a distance of 750pc (Matthews 1979). Classified as a Class 
0/ 1 borderline source, this young protostar is associated with a 
very energetic bipolar outflow (Codella et al. 2001). In addition, 
near-infrared images by Nisini et al. (2001) have revealed the 
presence of a collimated jet. Thus far, it is one of the best stud- 
ied Class 0/1 intermediate-mass sources. 



2. Observations 

We used the six-element IRAM Plateau de Bure Interferometer 
to observe IC 1396N simultaneously in three hyperfine transi- 
tions of CH 3 CN (7=13^12) at 239.064GHz, in the A- and 
E-state transitions of CH 3 OCHO (/=7->6) and (7=8^7) at 
90.188GHz, and in the continuum at 1.3mm and 3.3mm. 
Observations were made in configuration A (January 27, 2005), 
new A (February 12, 2006) and in a non-standard configuration 
(March 22, 2006) in conditions of excellent atmospheric seeing 
(0.3-0.5" FWHM) and good atmospheric transparency (pwv = 
1-4 mm). The spectral correlator was adjusted to observe the 
line transitions with a velocity resolution of 0.3 km s _1 and to 
cover the entire RF passbands (580 MHz) for highest contin- 
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Table 1. Millimeter flux densities, sizes, spectral indices and masses 



3.3mm (90.3 GHz) 


BIMA3 




BIMA2 








Cocoon 


41.86+11.9 


41.73+12.8 


41.73+14.3 


a (J2000) 


21 h 40 m 4Z84 


21"40 ra 41 s :86 


2lM0 ra 41 s :85 


21»40 m 41 ! :73 


21 tl 40 m 4f!72 


6 (J2000) 


58°16'01'.'4 


58°16'13'.'2 


58°16'11'.'9 


58°16'12'.'8 


58°16'14'.'3 


Size (") 


0.8" x 0.5" 


4.3" x 3.1" 


< 0.5" 


< 0.5" 


< 0.3" 


S (mJy) 


8 


16 


5.9 


1.5 


2.7 


1.3mm (237.6 GHz) 


BIMA3 




BIMA2 








Cocoon 


41.86+11.9 


41.73+12.8 


41.73+14.3 


a (J2000) 


21 h 40 m 4Z84 


21 h 4() m 41 s 85 


21 h 40 ra 41!86 


21 h 40 m 41 s :73 


21 h 40 m 41 s ;73 


6 (J2000) 


58°16'01'.'4 


58° 16' 13'.' 1 


58°16'11'.'9 


58° 16' 12'.' 8 


58°16'14'.'3 


Size (") 


0.8" x 0.4" 


4.5" x 3.1" 


0.4" x 0.2" 


< 0.3" 


< 0.2" 


S (mJy) 


30 


245 


35 


6 


10 




BIMA3 




BIMA2 








Cocoon 


41.86+11.9 


41.73+12.8 


41.73+14.3 


Mean Spec. Index 


1.4 






1.9 




Spec. Index (a) 


1.4 


2.8 


1.9 


1.5 


1.4 


Mass (M Q ) fl 


0.05 


0.4 


0.06 


0.01 


0.01 


/3 = a-2 


< 0.0 


0.8 


-0.1 


<0.0 


<0.0 



masses have to be scaled up by a factor of ~6 for a dust temperature of 20 K. 



uum sensitivity. The overall flux density scales for each epoch 
and for each frequency band were set by comparison to the 
known flux densities of MWC349 and 3C273 and estimated to 
be accurate to 10% at 3.3mm and 20% at 1.3mm. The resulting 
line and continuum point source sensitivities were estimated to 
8.4 mJy (2.5 MHz) and 0.43 mJy at 237.571 GHz, and to 2.5 mJy 
(2.5 MHz) and 0.12mJy at 90.250 GHz. The corresponding syn- 
thesized beams adopting uniform weighting were 0.45" x 0.37" 
at PA=-9° and 1.2" x 1.0" at PA=-11°. The data were cali- 
brated, mapped and analyzed in the GILDAS software package. 

3. Continuum data 

Figures la and lb show the continuum images obtained from the 
visibility data in the two frequency bands. The high sensitivity 
and the striking similarity of the many morphological features 
in the two figures, clearly testify to the presence of at least three 
bright continuum emission cores in the center of IC 1396N, all 
three very likely associated with the centrally peaked dust emis- 
sion source identified as BIMA 2 by Beltran et al. (2002). While 
the two weaker cores were not resolved by the interferometer, 
the primary core IRAM 2A is resolved in the 1 .3mm emission to 
an elliptical region of ~300 AUx 150 AU. This corresponds to 5x 
the size of the solar-system, and may be an accretion disk sur- 
rounding the powering source of the bipolar outflows mapped by 
Beltran et al. (2004) in various molecular transitions. The bulk 
of the mm-emission, however, is emerging from a much larger 
region centered on the triple system and oriented almost perpen- 
dicular to the molecular outflows (see Beltran et al. 2002). This 
oval-shaped region whose extent (~2800 AU) is very much the 
same at both observing frequencies, provides about 50% of the 
1 .3mm flux density detected at the JCMT (Saraceno et al. 1996). 

According to these results, we consider two different mod- 
els for the continuum emission: a large volume of cold dust (re- 
solved out by the PdBI) that surrounds (a) an envelope with sharp 
boundaries in which the three compact cores are embedded, 
(b) a confined region harboring multiple lower brightness cores 
from which we have detected the three most intense. Though 
the lack of sensitivity to large scale emission and emission dis- 



tributed over a large number of cores makes it difficult to argue 
against one or the other model, we favor the former, more sim- 
ple model of the dusty 'cocoon' at the border of which the three 
intense cores are embedded. We obtain excellent results by fit- 
ting such a model at 3.3mm and 1.3mm, which strengthens the 
'cocoon+cores' interpretation for the protostellar region associ- 
ated with BIMA 2. The fitted parameters of our model and the 
corresponding model image are shown in Table 1 and Figure 1 
(bottom). We cannot, however, totally discard the possibility of 
the existence of a relatively homogeneous cluster of lower mass 
cores surrounding the three bright cores detected. 

In addition to the intense clump BIMA 2, Beltran et al. 
(2002) detected two clumps in continuum emission, BIMA 1 
and BIMA 3. In spite of the much higher sensitivity of the ob- 
servations presented here, BIMA 1, the clump in the south-west 
of BIMA 2, was not detected in the 3.3mm continuum image. 
It is likely that the interferometer resolved out extended (< 5") 
low brightness continuum emission from this source. BIMA 3, 
the radio-bright source in the south-east of BIMA 2, is detected 
at the edge of our continuum images. The continuum emission 
towards BIMA 3 can be interpreted as compact emission emerg- 
ing from a region ~450 AU in size, embedded in a larger cloud 
of warm and dense dust and resolved out by the interferometer, 
accounting for about 70% of the continuum emission losses at 
1.3mm. 

3.1. Spectral Indices and Masses 

In Table 1 we present the primary beam corrected flux densities, 
sizes and spectral indices derived by fitting our model to the vis- 
ibilities in the uv-plane. The spectral indices range from 2.8 for 
the cocoon to ~ 1 .4-1 .9 for the compact cores. While the spectral 
index for the cocoon is consistent with optically thin dust emis- 
sion with p ~ 1, the value measured towards any of the compact 
cores is very low {a < 2). 

Different explanations may account for the low spectral in- 
dex measured towards the compact cores. One possibility is 
that the flux densities are a mix of dust and free-free emission. 
Beltran et al. (2002) reported VLA continuum flux densities 
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at 3.6cm and upper limits at 12cm for BIMA 2 and BIMA 3. 
Since BIMA 2 presents a complex structure at the high angu- 
lar resolution of our observations, we compare the flux densi- 
ties at centimeter and millimeter wavelengths in BIMA 3 only. 
We have extrapolated the emission at 3.6cm to millimeter wave- 
lengths assuming the spectral index of an isotropic ionized wind 
(S y oc v 6 ) and subtracted it from the PdBI flux densities. The 
resulting values were then fitted assuming optically thin dust 
emission. The ratio between the 1.3mm and 3.3mm emission 
can then only be fitted if we assume a spectral index of /? ~ 
for the dust opacity. One possibility is that the dust emission is 
optically thick, which is, however, not consistent with the mea- 
sured flux densities. With atj 3mm = 0.01crrT 2 g _1 , a gas and dust 
mass of ~1.1 M is needed to yield opacities >1 in BIMA 3. 
A moderately warm core (Td ust ~50 K) of 1.1 M would emit a 
flux density of 300 mJy at 1.3mm, i.e. a factor of 10 more than 
observed. Although our observations are not consistent with an 
optically thick core of ~450 AU, and despite a relatively low sur- 
face brightness temperature (~4K) at 1.3mm, we cannot discard 
the possibility that part of the core is optically thick. Clearly, an 
optically thick region would contribute to a lower value of the 
spectral index. 

Another possibility is to assume a different distribution of 
grain type in the compact cores and in the surrounding cocoon. 
Low values of ft (~0-0.5) are found in circumstellar disks around 
Herbig Ae and Be stars (Natta et al. 2004, 2007; Fuente et 
al. 2003, 2006) and have been interpreted as evidence of grain 
growth. According to this interpretation, the differences in the 
spectral index we measure in the cocoon and in the compact 
sources could hint at differences in the grain properties at dif- 
ferent spatial scales. These call for a scenario in which the com- 
pact components are circumstellar disks embedded in a dust- 
enshrouded enviroment. Similar differences in the spectral index 
between the envelope and the compact components of Class 
sources were reported by J0rgensen et al. (2007). 

To derive mass estimates for the different components (see 
Table 1), we have assumed a dust temperature of 100K for the 
cocoon and the cores. Although the dust grain properties are 
very likely to be different, we have based our mass estimates 
on a uniform dust emissivity index K\ ^ mm = 0.01cm _2 g _1 . All 
the compact cores have masses much lower than one solar mass, 
possibly suggesting that the PdBI has resolved out a large frac- 
tion of the envelope. We can compare the mass of the cores with 
those derived towards hot corinos and intermediate mass hot 
cores. After correcting for the distance, the 1.3mm flux density 
towards IRAM2A is a factor of 3-4 larger than the one mea- 
sured towards the hot corinos IRAS 16293-2422 A/B (Bottinelli 
et al. 2004a) and testifies to the presence of a more massive star. 
However, IRAM2A is lOx weaker than the intermediate mass 
hot core NGC 7129-FIRS 2 suggesting that IRAM2A is less 
massive, although both protostars are located at about the same 
distance. Despite the fact that IC 1396N has a luminosity simi- 
lar to NGC 7129-FIRS 2, the envelopes surrounding these stars 
very likely correspond to very different star formation regimes. 
While NGC 7129-FIRS 2 at prima facie hosts a single massive 
core (pending higher resolution observations), IC 1396N defi- 
nitely hosts a cluster of lower mass stars. We argue that the most 
massive component of the IC 1396N cluster is very likely the 
precursor of a Herbig Ae star. 

4. Molecular line data: CH 3 CN 

CH3CN (1 3^t — >\2k) line emission has only been detected to- 
wards IRAM2A, the most massive core in the IC 1 396 N proto- 



cluster. In addition to the CH3CN lines, several CH3C2H lines 
appear in the spectrum testifying to the presence of a core, rich in 
complex molecules (see Fig. 2). The integrated intensity map of 
the CH3CN (13jt— >12t) transition shows that the emission arises 
in a region of radius ~0.8" (600 AU) around the massive hot core 
IRAM 2A and with a good positional coincidence between line 
and 1.3mm continuum emission. 

Comparing the integrated emission with the peak flux den- 
sity, we estimate that ~50% of the emission is coming from a 
compact region (radius <150 AU) around the hot core while the 
remaining 50% arises from a more extended and cooler com- 
ponent in the immediate vicinity. Bottinelli et al. (2007) report 
about a similar case in NGC 1333 IRAS 4A. Although the at- 
tained sensitivity is not deep enough to draw firm conclusions on 
the small-scale spatial distribution of the CH3CN emission, the 
shape of the emission suggests that CH3CN is tracing the signa- 
ture of gas and dust in the hot core and possibly arises from a nar- 
row boundary layer at the outflow-core interface. The marginal 
or non-detection of CH3CN towards the weaker cores is possibly 
subject to the same sensitivity limitation. 

We have tried to fit the rotational diagram of the CH3CN 
(13^ — > 12^) transition towards the emission peak. The K-ladder 
cannot be fitted with a single rotational temperature suggest- 
ing that the CH3CN emission is optically thick in the hot core. 
Also, we cannot exclude that some K-components are blended 
with other lines. Given all the uncertainties, we derive a CH3CN 
abundance using the temperature inferred from the rotational di- 
agram. We obtain a rotational temperature of 240±150 K and 
Af(CH3CN)= 7.1 10 14 cirT 2 towards the peak, which implies a 
CH3CN abundance of 2 10~ 10 . This abundance is similar to the 
ones determined by Bottinelli et al. (2004b) in the hot corinos 
NGC 1333 IRAS 4 A and IRAS 16293-2422 but a factor of 10 
lower than the one measured by Fuente et al. (2005) towards 
NGC 7129-FIRS 2. This low abundance of CH3CN also argues 
in favor of IC 1396 N being a cluster of low /intermediate mass 
stars. 

5. Summary 

We report on interferometric 3mm and 1.3mm continuum and 
CH3CN (13t— >\2k) line observations of the intermediate-mass 
protostar IC 1396N in the most extended configuration of the 
PdBI. Our observations provide a spatial resolution of ~250 AU 
at the distance of IC 1396 N and reveal the existence of a cluster 
of cores in this IM protostar with unprecedented detail and qual- 
ity. The cluster in the interior of IC 1396N comprises at least 
three cores, the very likely precursors of two low-mass and an 
IM star. Because of sensitivity limitations, we cannot exclude 
the presence of further low brightness cores in the BIMA 2 re- 
gion. Our continuum observations also show that the spectral 
index at millimeter wavelenghts changes with the spatial scales 
traced by the interferometer. The variations in the spectral in- 
dex are very likely driven by grain properties changing with gas 
density and by the existence of compact regions in the hot cores 
with optically thick emission. 

CH3CN (13— >12) line emission has only been detected to- 
wards the most massive core IRAM 2A, but our sensitivity limit 
is not low enough to infer a different abundance in the other 
cores. The CH3CN emission is slightly more extended than 
the associated dust emission and possibly arises from a region 
heated up by the molecular outflow. 
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Fig.l. Continuum map of IC 1396N at 90.250 GHz (top) and 
237.570GHz (middle). The lcr noise levels are 0.12mJy/beam 
(0.02 K) and 0.43mJy/beam (0.06 K), respectively for the 
3.3mm and 1.3mm continuum. The ellipse in the lower left cor- 
ners shows the synthesized beam, the dashed ellipse in the mid- 
dle panels shows the FOV of the 1.3mm receivers. The cen- 
ter of the dashed ellipse marks the position of the phase center 
at a = 21 h 40 m 41^71 and 5 = 58°16'12'.'8 (J2000). The figure 
(bottom) shows the 1.3mm reconstructed and primary beam cor- 
rected model image. BIMA 1 (open star) is to the S W of BIMA 
2. 
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Fig. 2. CH 3 CN (13* -» I2 k ,k = 0...5) and CH 3 CCH (14* -> 
I3 k ,k = 0. . .4) lines at 239.064 GHz observed with the PdBI 
towards the core 41.86+11.9 (= IRAM2A). The intensity scale 
is in units of mJy and the spectral resolution is 2.5 MHz (= 
3.1 kms -1 ). The one sigma noise is 8.4 mJy. 
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Fig. 3. Continuum emission at 1.3mm superposed on the map 
of the integrated intensity flux of the k = 0,1,2 components of 
the CH3CN 13*— >12* line (contours). Levels for the continuum 
emission are from 3<x to 9cr in steps of l<x = 0.12Jy km s _1 . 
Significant CH3CN emission is detected towards the most mas- 
sive core only. 



